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ABSTRACT 

We discuss the correlation function of hadronic currents in the Schwinger model at 
finite temperature T. We explicitly construct the retarded correlator in real time and 
obtain analytical results for the Euclidean correlator on a torus. Both constructions lead 
to the same finite temperature spectral function. The spatial screening lengths in the 
mesonic channels are related to the dynamical meson mass m = e/y^ and not 27rT even 
in the infinite temperature limit. The relevance of our results for the finite temperature 
problem in four dimensions is discussed. 
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1. Introduction 

With the possibility to use ultrarelativistic heavy ion coUiders to probe hadronic matter 
at very high temperatures and/or densities it becomes imperative to get an improved 
understanding of finite temperature gauge field theories. A variety of approaches have 
been devised ranging from first principle lattice simulations to approximate estimates 
in simplified models 0. 

Just as in conventional plasma physics, a central problem in finite temperature QCD 
is to find the excitation spectrum as a function of temperature, density, external fields, 
etc. . The spectral information in some particular channel a is contained in the spectral 
function (Tq. Its knowledge would, at least in principle, allow us not only to calculate 
responses to external probes, but also production rates, multiplicities, etc. . 

While the spectral function at zero temperature contains poles corresponding to bound 
states and resonances and cuts corresponding to particle production thresholds, its inter- 
pretation at finite temperature is not always clear. It is usually assumed that the spectral 
function at finite temperature is dominated by hadronic quasiparticles in the low tempera- 
ture phase and various collective excitations in the high temperature phase. However, the 
nature and quantum numbers of these collective excitations are not very well understood. 

To determine the spectral function in a theory like QCD from first principles remains 
an unsolved problem. Besides the purely technical/numerical problems, the present lattice 
simulations of finite temperature systems can primarily provide information about space- 
like correlators only. The spectral function at finite temperature, which is related to 
time-like correlators, can only be obtained if one makes very specific, and unprovable, 
assumptions about its analytic structure. 

Lacking suitable lattice techniques, other methods have been used to get spectral in- 
formation. Examples are finite temperature QCD sum rule calculations and phenomeno- 
logical models (usually inspired by models for the zero temperature vacuum), but here one 
is forced to assume a special shape for the spectral function (shifted poles and cuts) and 
then proceed to estimate some unknown parameters like resonance masses and thresh- 
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olds (see for a review). Another approach advocated by two of us 0], is to use the 
method of dimensional reduction to calculate certain static correlation functions at finite 
temperature. 

The questions we would like to raise in this paper are twofold: 1. Is it meaningful to 
talk about spectral functions in interacting field theories at finite temperature? 2. To what 
extent can we learn about them from finite temperature Euclidean correlators? We will 
not attempt to answer these questions for QCD, but use a model theory simple enough 
to obtain answers. 

We consider the Schwinger model at finite temperature. In section 2 we explicitly 
construct the real time correlation function for pseudoscalar and scalar sources (analogue 
of e^e~ production channel) at zero temperature and derive an explicit relation for the zero 
temperature spectral function. In section 3, we extend the analysis to finite temperature, 
by considering the retarded response function, and derive an explicit expression for the 
finite temperature spectral function for the Schwinger problem, using operator techniques. 
In section 4, we evaluate the scalar and pseudoscalar correlator for the Schwinger model 
on the torus. Exact analytical results are obtained with the help of the properties of 
the Jacobi theta-functions. The results are in full agreement with the naive bosonisation 
approach. Issues related to the temporal and spatial screening masses are discussed in 
section 5. In section 6, we show that the analytical continuation of the full nonperturbative 
Euclidean correlators yields temperature dependent spectral functions that are identical 
to the ones derived from the real time approach. Our conclusions are summarized in 
section 7. Helpful calculational details are given in two appendices. 

2. Feynman Correlators 

The Schwinger model is quantum electrodynamics [^] in one space and one time di- 
mension defined by the Lagrangian density 

C = \Fl-i,tpil^, (2.1) 
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where F^jy = d^A^—d^A^j, is the Abehan field strength and p= •y^^d^—ieA^) the covariant 
derivative. It is well known that at zero temperature this model is equivalent to a theory 
with a single massive and free meson with mass m = e/ a/tt where e is the electric charge 
of the electron [^]. The meson field is related to the polarization density induced by the 
electrons. 

The mesonic Green's functions in the vacuum can be constructed using the bosonisa- 
tion techniques discussed by Coleman [0. In this paper we consider the following scalar 
and pseudoscalar correlators, 

Pf{x) = (T(^^75^(a;)^275^(0))), (2.2) 

which can be readily evaluated by using the operator identity 

± 75)z^(x) = - : e^^^v^*(^) :, (2.3) 

where is a free massive meson and C = 0.577 is the Euler-Mascheroni constant. Normal 
ordering with respect to the vacuum is denoted by : :. The field (pix) can be decomposed 
into positive and negative frequency components as 

0(x) = / (e-*^-^°+^^^'ap + e+'''-^°-'P^"4) = 0-(x) + 0+(x), (2.4) 

with Ep = y/p'^ + m^. The product of normal ordered operators obtained by inserting 
( p.3|) in (|2.2| ) can be evaluated with the help of the following relation {exfi = ±1) 

which can be established using the Baker- Campbell-Hausdorff formula. Inserting (|2.3|) 
and ( |2.4| ) into (|2.2| ) and using repeatedly ( p.5| ) yield 



Pp[x) = {ipijfsmh{-4mAF{x,m^)), (2.6) 
where the Feynman propagator is given by 

AF{x,m') = e{x')[r{x),<p^m-0i-^'M^i^),rm, (2.7) 



and the fermion condensate is 

m ^ - (^) . (2.8) 

Similarly, 

Sp{x) = (^V^)2cosli(-47riAi7(x,m2)). (2.9) 
These results were first obtained by Casher, Kogut and Susskind By expanding the 



hyperbolic functions in (|2.6|) and (|2.9| ) we observe that the pseudoscalar correlator is a 
sum of exchanges of an odd number of massive mesons while the scalar correlator ( p.9[ ) is 
given by exchange of an even number of mesons. The scalar field carries odd parity. In 
addition, the scalar correlator is nonvanishing and equal to {'ipip)'^ in the x oo limit as 
expected from the cluster decomposition. The nonvanishing of {ipipy in this model was 
first discussed by Loewenstein and Swieca 

The spectral function at zero temperature follows from (^.21) by inserting a complete 



set of physical states and making use of translational invariance, 

POO 

Pp[x) = -i dm^ap{m^)AFix,m^), (2.10) 
Jo 

where we have defined (g° > and = m^) 

apiq') = 2nY^5'iq-pn)m'^tl,m\n)\'. (2.11) 

n 

A comparison of (|2.6|) with ( p.lO|) yields the following spectral function for the pseu- 
doscalar channel 

oo 22fc+l 2k+l I. 2fc+l 



i=i i=i 



(2.12) 



which was also obtained in [|T^. The first term corresponds to the meson pole with a 
strength determined by the quark condensate. The sum is over the 3, 5, ■ ■ ■ , oo-particle 



phase space with a couphng strength proportional to the pole strength. The spectral 



function ( 2.12 ) is just the sum over the cuts of the Feynman diagrams in the expansion 
of the correlator Pp in powers of A^(a;). A similar expression can be derived for Sf{x) 
involving an even exchange of pseudoscalar mesons. 



3. Retarded Correlators 



The finite temperature spectral function can be obtained from the retarded correlator. 
For that consider for instance the pseudoscalar correlator 



Pn{x) = ^(x°)([^^75^(x), V'^75^(0)])^, 



(3.1) 



where {■ ■ ■)i3 denotes the Gibbs average with respect to a heat bath of mesons of mass 



m = e/ -\/7r. The expectation value in (3J.) can be calculated using the operator identities 
( ^73|) and (|2.5|) . In a heat bath the expectation value (0+0^ )/3 does not vanish so that 



(exp (iie^0+(a;) + zeo0+(O)) exp {ie^cj) {x) + zeo0 (0)))/3 = 

exp (-eoeo(0+(O)0"(O))^ - e.M'P^{^)'P'{^))p-^oe.{<P^{m-{x))p - e,eo(0+(x)r (0))^ 

(3.2) 

With this in mind, the retarded pseudoscalar correlator becomes 
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47r2 



sinh (47r([0-(x), 0+(O)] + 47r(0+(x)0"(O))^ + 47r(0 
sinh (47r([0-(O), 0+(x)] + 47r(0-(x)0+(O))/3 + 47r(0^ 



-(O)0-(x));3 



(3.3) 



In the heat bath 



[0-(x),0+(O)] + (0+(x)0-(O))^ + (0+(O)r (a:))/3 



(3.4) 



where np = l/{e^^p — 1) is the Bose- Einstein distributions for the massive bosons. Using 
3) in (O) gives 



sinh (^47r J ^e-'P--6{p'-m'){e{p^)+nS 
sinh Utt J ^e-'P-^S{p' - m2)(^(-/) + n^] 



where 



,(3.5) 



(^•W. - -—e-"~ . -_ exp ^- y j (3.6) 



is the temperature dependent chiral condensate, as calculated in [jTT], [T^, |T3|. A similar 
result holds for the retarded scalar correlator S'ij(x) with the substitution sinh — cosh in 
( p.SI ). Micro-causality can be exhibited explicitly in ( p.5[ ) by rewriting the result in the 
form 

(fp 



Pr{x) = sinh(27rA«(x,m2)) cosh ^(1 + 2^^) ImA^(p, m^) e^f-j , (3.7) 

where Aij(x, m^) is the retarded propagator 

A^(x,m2) = ([0+(x),0-(O)] + [0-(x),0+(O)]) , (3.8) 

and /S.p{p,im?) the Fourier transform of the Feynman propagator ( p.7| ). Again, the re- 
tarded pseudoscalar propagator in the heat bath corresponds to the exchange of an odd 
number of pseudoscalar mesons with temperature dependent strengths. The temperature 
dependent strengths follow from the bubble insertions at the vertices. Fig. la shows the 
temperature insertions on the vertices, and Fig. lb the diagrammatic expansion of ( p.7|) . 
A similar result can be derived for the scalar correlator. The answer is 

Sr{x) = sinh(27rA^(x,m2)) sinh ^(1 + 2^^) ImA,.(p, m^) e'^'^^ (3.9) 

and involves an exchange of an even number of pseudoscalar mesons, as expected from 
parity. Fig. Ic shows the corresponding diagrammatic expansion. 

Using translational invariance, the retarded correlators can be related to the pertinent 
finite temperature spectral function. In the pseudoscalar case we have 

Pr{x) = J d\ ^e^<?Vp(g°, q') (3.10) 



where the temperature dependent spectral function reads 



A,B 



(3.11) 



The free energy of the system is denoted by F and q\ = m?. The spectral function for 
the retarded correlator is obtained by expanding the sinh in (|3.5| ) in a power series. To 
proceed the following identity is useful 

2fc+l 2fc+l 2fc+l 

n + - n + = 2sinh(/3//2) n e^'^-'/'n,., (3-12) 

i=l i=l 1=1 

where p° = J2iP^- The proof of this identity follows immediately if one realizes that 



n (1 + n - n n d + = , — • (3.13) 

By comparing the expansion of the correlator and (|3.1CI|) we can read the spectral function 



00 r)2k+l 2k+l 2k+l 



(3.14) 



A rerun of the above argument for the retarded scalar correlator gives (following the 
substitution sinh — > cosh in (P3|)) 



as{q\q^) = 47r(V'V')| sinh(/5go/2) 

00 o2fc 2k „ 2fc 

X E Tf^ n / - m^) e^l-^'l/^n,^5^(g - Ep.)- 



(3.15) 



The spectral functions at finite temperature are just the zero temperature spectral func- 
tions with temperature dependent strengths and thermally weighted phase space. The 
thresholds remain temperature independent. This result was expected, since the bosonised 
version of the model is a model with a free single massive meson. 
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The above results are consistent with the conventional cutting rules in real time. 
Specifically, we have for the scalar {pi = — p\) channel 



a5(g°,g^)=47r(^^)^ 



2Ei2E2 



X 



- El - E2)((l + r2i)(l + n2) - nma) + + + E2){nin2 - (1 + ni)(l + ^2)) + 
- El + E2)(n2(l + ni) - ni{l + n2)) + + Ei - E2)(ni(l + ^2) - n2(l + rii))) + 



(3 



where the dots refer to the processes with 4, 6, ■ • ■ pseudoscalar meson exchanges. For the 
pseudoscalar channel the first two terms in the expansion of the sinh-functions in (3.5) 
result in (pg = ~ p\ — pi) 



— 1 f)7r — 

ap{q\ q') = 47r(^V)j5(g' - m') + -^ml 



dp\dpl 



3 '^^'PJ 2E12E22E3 
-E,-E2- E,) - + E, + E2 + E^Mil + ni)(l + n2)(l + ^3) - 71^2713) 

- El - E2 + E3)((l + ni)(l + 722)^3 - ^1722(1 + ^3)) + (3 ^ 1) + (3 ^ 2) 

+ El - E2 + E3)(ni(l + n2)n3 - (1 + ^1)712(1 + %)) + (2 ^ 3) + (2 ^ 1)) + ■ ■ ■ , 

(3.17) 

where the dots refer to the processes with 5, 7, • ■ ■ pseudoscalar meson exchanges. Note 
that the pole contribution in the pseudoscalar channel remains temperature independent, 
only its strength is affected by temperature. Fig. 2a shows the various cuts contributing 
to (|3.16| ) and Fig. 2b the various cuts contributing to ( |3.17|) . The thermal contributions 
to ( |3.16| - |3rT7D reflect on the exchanges of the sources with the heat bath. Incoming arrows 
are weighted by n, outgoing arrows are weighted by (n + 1). Lines on the left hand side 
of the vertex denote absorption into the heat bath, lines on the right hand side of the 
vertex denote emission from the heat bath. The lowest contribution to the scalar channel 
is consistent with the result discussed by Weldon in four dimensions for scalar 0^ 
theories. The next to leading order contribution to the pseudoscalar channel is consistent 



with the result discussed by Fujimoto et al. |15] in the context of scalar 0^ theories. The 
results ( p. 141) and (|3.15| ) summarize all the cuts for scalar 0" theories with n = 3, 4, ■ ■ ■ , cxo. 
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Their extension to four dimensions is immediate. The singularity structure of the spectral 
function at finite temperature is much more involved than the zero temperature one. Its 
structure is physically transparent when expressed in the energy variable g° (in the heat 
bath frame) as opposed to the invariant variable g^. 



4. Euclidean Correlators 

Most calculations in finite temperature field theory are performed using the imaginary 
time or Euclidean formulation. At the perturbative level it has been well established how 
these results are related to similar calculations in real time. However, in the previous 
section we obtained the full nonperturbative correlation function based on operator tech- 
niques in real time and it is a nontrivial question whether it can be reproduced by an 
Euclidean calculation. In Euclidean space, finite temperature corresponds to imaginary 
time with a periodicity f3 = 1/T. The infrared divergences in the Schwinger model are 
regulated by putting the space axis on a circle with circumference L that at the end of 
the calculation is taken to infinity [0. Thus, we are led to study the Schwinger model 



on a torus. 

We now outline the calculation of the scalar and pseudoscalar correlation functions on 



a torus. We will follow the methods developed in ||16| , and and use the notation of 
the latter reference when needed. In the following, all expectation values are defined on 
the torus except when stated otherwise. The starting point is the generating functional 
for (disconnected) fermionic Greens functions, 

k 

Z[q,f]] = (n(r7|7/^p)(V^p|r7)e-/'^'"'^'^^(")^(""J')''(j')det'(i^))A. (4.1) 
p=i 

The fermionic Green's function S{x — y), with proper boundary conditions, is derived in 
Appendix A, and the gauge- field average {■ ■ ■)a is with respect to the usual Maxwell action 
— ^E"^. The functions ipp are the zero modes, and as explained in [|ll|, in the topological 
sector corresponding to a non zero electric flux $ = 27ik/e, there are precisely \k\ zero 
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modes, all with cliirality $/|$|. In this sector the gauge potential can be decomposed as 

^ = jho + di(l) + do\-^x\ (4.2) 
Ai = ^/ii -9o0 + 9iA, (4.3) 



where the are constant and V = (3L is the area of the torus. The explicit expression 
for the p^^ zero mode (p = 1, 2, ■ ■ ■ |A;|) with chirality a = ± can be found in ||11||: 



-aip{x) 



g27ri[/iixV/3-fca;"xV^] Q 



xyL + {p-i-hi)/k 

kx^/(3 + ho 



{0,i\k\T) (4.4) 



where is the generalized theta function as defined in |18|. The determinant in (|4.1|) is 
over the states orthogonal to the zero modes and is given by |TT | 



det(i^) 
for k = and by 
det' (iJl)) - 



T]{iT 



2\k\ 



-0 



1 + ho 



(0,^r) 



/ 2 

exp ( j (f' X (j){x)^(j){x) 



\k\/2 



exp j £x<p{x)U<p{x) + J 



X0[X) 



(4.5) 



(4.6) 



for k ^ 0. Here r]{iT) is Dedekind's eta function. 

As in the previous section we study the connected scalar and pseudoscalar correlation 
functions defined byQ 

Se{x) = {^tlj{x)^^{0)), 
Pe{x) = (V^^75V'(a;)?^75V'(0)), (4.7) 

where the average is according to the partition function (4.1). It is convenient to use a 
chiral basis for the spinors and the fermionic Greens functions defined by 



^,{x) = Peip{x) = -(1 + e75)V'(a;), 



(4.8) 



where again e = ±1 and similarly for S^e'{x — y)- The above correlation functions can be 
expressed as 



Se{x) 
Pe{x) 



2C+_(x) + 2C. 



2C+_(x)-2C++, 



(4.9) 



^Our conventions for the 7-matrices in Euclidean space are : [7^, 7^]- 



1^1., it = 7m and 75 = 17172. 
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where the Cee'{x) are defined by 

C«,(x) = {^PMx)^Pe'^iO)). (4.10) 

and we have simphfied ( ^4.7| ) by using the relations 

= = C++(x). (4.11) 

From ( [4.ip and the chirahty of the zero modes it follows that in fact only k = and 
k = ±2 contribute to the correlators (|4.7| ). In the trivial sector only C+^{x) and C_+(x) 

are non- vanishing whereas C++(x) and C (x) are non- vanishing in the sectors k = 2 and 

/c = —2, respectively. 

We start with the trivial sector where the relevant non-vanishing correlation function 
reads 

C+_(x) = (^P+V^(x)V^P_z^(0)) = -^(det(2^) S^4x)S+4-x))a, (4.12) 

where Z = Z[0,0] is the partition function on the torus. To evaluate the functional 
integral over the gauge field we again follow []TT| and write 

(4.13) 

where 0m and Am are the Fourier components of the fields and A subject to the condition 
/ d'^x(f){x) = J d'^x X{x) = 0. Furthermore, the Jacobian J does not depend on the 
dynamical fields and cancels between numerator and denominator in ( [4.12| ), and the 
integration over the pure gauge A is trivial. Substituting the expression (A. 11) for the 
fermionic Greens function, and the explicit Fourier expansion of the generalized theta 
function, it is now straightforward to perform the /i-integrations to get 



n {-')^ 1 



^'i(O) 



^i(^) 



2 



^g-2<^(.)+20(O)^^^ (4.14) 



where z = (x° + ix^)/P, and (■ ■ ■)s denotes the connected Greens functions with respect 
to the measure HmT^o '^0m exp(— 6*5) induced by the bosonic action 

Sb = \ I d^x(t){x)U{U - m^)0(x). (4.15) 



2 
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For simplicity of notation we have introduced the Jacobi theta function 'i9i(a;) = —6ii{x, ir). 
Finally using the identity ||19| 



??;(0) = 2TXT]{iT) 



(4.16) 



we obtain 



1 



{e 



-20{a;)+2(/-{O) 



The integral over is Gaussian and can be expressed immediately as 

where the two-point function K{x, y) is defined by 

K{x,y) = {x 



□ (□ - m 
Here Am{x) is the Greens function 

1 



^ \y) = ^[Am{x - y) ~ A{x - y)]. 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



and A its massless counterpart. The full massive Greens function is obtained by adding 
the contribution from the constant mode, i.e. 



Ara{x -y) = Arn{x - y) - 



Vrn? 



(4.21) 



Thus we have 



^g-2(/.(a:)+2</.(0)^^ ^ ^ie'^K{x,x)-iTTl\,^{x)+4:-KA{x)-^^ 



(4.22) 



Note that ( [4.17|) apparently has a free massless singularity at = since "dii^z) is odd. 



However, the expression (4.22) for the massless propagator was studied by Kronecker [|19 
with the following result 



A(x) = — log 



—e V 



(4.23) 
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so this singularity cancels^. 

Combining equations (|4.17|) , (|4.22|) and ( [4.23|) and using the result 



(^^)^ = _|e-^|r7(^r)|2e2^'^("'") (4.24) 



from 11 1 we finally get 



C+_(x) = ^|^(,^)|4e4e2X(x,x)-4.A™(x) ^ 1 ^^^^ 2 g4.A„ (x) _ ^4_25) 

Note that this contribution to the scalar correlator does not saturate the cluster decom- 
position according to which 

jim (^?/'(a;)V^?/'(0)) = l^^i^f. (4.26) 

From this it is clear that the missing piece must come from the k = ±2 sectors which we 
now evaluate. 

According to the partition function ( [4.1|) , the correlation function C++(x) reads 



C++(x) = I < det'{tI^)\Mx)MO) - Mx)MO)\' >A . (4.27) 
The average over and hi is performed exactly in Appendix B resulting in 

o 7r 

C^4^) = ^l^(^r)re-^-^|^i(z)r(e^^(^H^^(°))5, (4.28) 



where we have used that the partition function is equal to Z = l/(|r7(zr)pv^2r). The 
integral over can be carried out in exactly the same way as in the k = sector. The 
result is 

(^^2<p{x)+2<l,{0)^ ^ ^ g4e2K(a;,a;)+47rA„(x)-47rA(x)+-^_ ^^^29) 

Again we use the Kronecker-identity for the massless propagator which leads to a cancel- 
lation of the ^—functions. Finally 

C++(x) = i_|,^(^r)|4e^^'^("'")+^"^'"(") = ^(V^7/^)Je-^"^'"("). (4.30) 
P ^ 



^The cancellation of the singularity is expected because this model does not have massless excitations. 
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Combining (^4.25|) and ( [4. 301) we get the remarkably simple results 



Se{x) = (?/'^)^cosli[47rAm(x)], 

Pe{x) = {^^)lsmh[AnA^{x)], (4.31) 

where Se and Pe are the Euclidian counterparts of ( p.9| ) and (|2.6|), respectively. These 
equations are derived for a finite circle, but as shown in the L ^ oo limit of {ipip)f3 is 
well defined and exactly equal to the temperature dependent condensate (ipip) p appearing 
in section 2. These expressions were expected from the bosonisation approach described 
in the preceding sections. Indeed, ( [1.31| ) are the natural versions of the free massive boson 
correlators on the torus. 



5. Screening Lengths 

Lattice simulations in four dimensions in QCD have focused on the Euclidean corre- 
lation functions which particular emphasis on screening lengths I^O], It is interesting 
to note that the spatial screening lengths following from ( [4.31j ) are controlled by the me- 



son mass m = e/y^ and not 27rT whatever the temperature. Indeed, for large spatial 



separation and at low temperatures (T <^ m/27r) ( [4.31| ) reduces to 



PEiO,x') ~ 2{i^ij)l Ko{m\x'\), (5.1) 

where Kq is the McDonald function. We have used the approximation 

1 r fin p^P^^ ( 2 \ 1 

valid for T <^ m/27r, and dropped terms of order e"™'/"^. In ( |5.1| ) the first result reflects 

on a cut corresponding to the exchange of two massive pseudoscalarsQ, and the second on 

■^Note that the pre-exponential factor is essential in determining the correct singularity structure of 
the spectral function. We believe this to be the case also in four dimensions and at low temperatures 
with 2; "'"I) ^ -ftri(TO|xi|). 
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a pole. For temperatures such that T ^ m/27i, {\i.31f) reduce to 



'ttT 



m / 



Pe{0,x') ~ 2{^^p)l _ e-'-l-^ (5.3) 



^ \ m 

where we have used Am(0,x^) ~ — Te~"^'^^'/2m instead of ( |5.2| ). We conclude that for 
any temperature, the screening length in the pseudoscalar channel is m and in the scalar 
channel 2m, whatever T. What is interesting to note is that for T ^ m/27r the finite tem- 
perature Euclidean correlator along the spatial direction, no longer correctly reproduces 
the singularity structure of the real time spectral function. 

If we were to dimensionally reduce the Schwinger model, we would expect the static 
correlators such as the ones described here to display screening lengths of the order of 27rT, 
reflecting on the free field contribution to the correlators at high temperature. Indeed, 
for a free spatial massless fermion {un = {2n + 1)7tT), 

5(0..') = -iE/|^^-'"'^^^^ = -g^fST. (5.4) 

p ^ J 2-71 uj^+p^ 2p smh(^) 

after contour integration. At large distances and/or large temperatures, only the lowest 
Matsubara modes ivrT contribute to (|5.4|) , resulting in 



S{0,x^) ~ -iTsgn(a;i)7^e-"^l^'l, (5.5) 

in agreement with the leading behaviour from dimensional reduction. The free fermion 
contribution to the scalar and pseudoscalar Euclidean correlators are 



in disagreement with ( |5.3| ). Why is that? The reason is due to the fact that the dimen- 
sionally reduced version is plagued with infrared divergences that makes the result (|5.6|) 
inapplicable. The model offers a counterexample to by now naive intuition. 
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The behaviour of ( [4.31| ) along the temporal direction at low temperature T -C m/27r, 
follows from the following form of Am{x^, 0) 



A^(x°,0)= i e +e e e 



The singularity structure in the scalar and pseudoscalar correlator following from (|5.7|) 
at very low temperature agrees with the one derived from the spatial asymptotics, as 
expected from 0(2) invariance in the zero temperature limit. The temperature corrections 
mostly affect the strengths of the singularities. These corrections are suppressed by e~'^l'^ . 

In the temperature range T ^ m/2Ti^ the behavior of ( |4.31| ) along the temporal 
direction is given by 

Se{x\ 0) ~ cosh f ^ - ln(2 - 2cos(^x°))) ~ 4 ^ ^ 



...... I - ^-^v^- ,,j 



P^(a:",0) ~ -(^^)^sinh — - ln(2 - 2cos(— x")) ~ -— — 



/3 

(5.8) 



We have used the fact that at high temperature [|TT 



{iPiP)p ~ -2re-'^. (5.9) 

The exponential growth in the hyperbolic functions in ( ^.8|) is balanced by the expo- 
nential fall off in the condensate. To compare this result with the free massless fermion 
contribution at finite temperature, we recall that for a free temporal massless fermion 
{Un = {2n + l)7rT) 

f3\^J 2n' ul+p^' - 2/5 sin(^)' ^^-^^^ 

where the last relation follows by contour integration with x*^ 7^ modulo j3. This 
immediately shows that at high temperature the scalar and pseudoscalar temporal cor- 
relators (|5.8| ) reduce to the free fermion behavior^ By squeezing the temporal direction 



4 



Note that the free results (|]J) and ( 5.10| ) in their exact form are analytical continuation of each 



other, even though the boundary conditions are not interchangeable for fermions. 
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in Euclidean space (high temperature) the correlator in the temporal direction picks up 
exclusively the free field part. This is not the case along the spatial direction. 
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6. Analytical Continuation 



We now show that the Euclidean correlators ( |4.31| ) give the same spectral functions as 



the ones derived from the retarded correlators. The derivation is fully nonperturbative. 
The Euclidean correlators at finite temperature are periodic. This can be readily seen by 
writing ([4.7|) in a Hamiltonian formalism. For instance {x = (r, x^)) 

PE{r,x') = Tr (e-^(^-^)(0(r) ^ij,tP{x) ^^75^(0) + ^(-r) ^i^^^O) i^^l^H^))) ■ (6-1) 

Using ^p{T + P) = —e'^^ il){T)e~^^ , it follows immediately that Pe{t + 13-, x^) = Pe{t, x^). 
Similarly for the scalar correlator. Thus Pe{t,x^) has a discrete Fourier representation 
{un = 27mT) 

PE{r,x') = ^^e-^-"^PijK,xi). (6.2) 

P n 



For < r < /9, the Euclidean correlator ( |6.1| ) can be rewritten in the form 

PE{r,x') = I dp^P'^'^Y.^ip' - (Pb -pi))e-^(^-^-^)e-(^--^-)-| < A|^^75^(0)|i? > |^ 



A,B 



(6.3) 



following the insertion of a complete set of physical states, = m^. Inserting (|6.3|) and 
the identity in frequency space into ( |6.2D , give 



PE{r,x') = 4e /^e— (6.4) 

The sum over the Matsubara frequencies lOn can be carried out explicitly, using the relation 
(for the pinch see below) 

-Ve"''^"^ = = i —— - f6 5) 

/3 V Jc2Tciz-LU e*^^ - 1 ■ ^ ■ ^ 



The contour C encloses the real axis counterclockwise. Inserting (|6.5|) into ( |6.4|) and 
unwinding the contour with < r < /3 in mind give 
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We have deliberately disregarded the pinch singularity in (|6.5| ) at = 0, as it gives zero 
contribution in ( |6.6|) since ap{0,p^) = 0. 

For L — > cxD the bosonic propagator Am{x) can be written as 

^m{r, x') = ~J: P ' . (6.7) 



The sum over n can be converted into an integral as in (|6T 

1 f dz f dp qW^^^ g«2-r 

Again, the contour C encloses the real axis counterclockwise. C can be unwound to give 
contributions only from the poles ±iy/p''^ + on the imaginary axis in the 2;-plane (with 
no pinch in this case). Thus 



The lack of manifest periodicity in r in ( |6.9[ ) is due to the fact that while unwinding the 
contour C above, we have explicitly restricted r in the range [0, (3] to allow convergence 
on the circles at infinity. The spectral function follows from (|4.31|) by expanding the 
sinh-function in powers of A^, using ( |6.9|) and comparing with (|6.6| ). The result is in 
complete agreement with ( p.l4| ) as obtained from the retarded correlator. This shows 
that the full real time spectral function can be reconstructed from the nonperturbative 
Euclidean correlator at finite temperature. 

In general, the full analytical knowledge of the finite temperature Euclidean correla- 
tor is a " dream-stufF' in interacting quantum field theories in four dimensions. We are 
therefore left with approximate, and usually numerical, estimates of the these correlators. 
Given this fact of life, how do we go about reconstructing the true spectral function? 
More specifically, given some knowledge of the l.h.s. of ( |6.6| ), how do we go about the 
r.h.s. of (|6.6D? 

At zero temperature, the conventional wisdom is to look at the asymptotics of the 
Euclidean correlators, under the assumption that the leading singularity in the spectral 
function is well separated from possible background. At low temperature T ^ m/27r, 
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this rationale seems to hold for the spatial asymptotics while most of the zero tempera- 
ture structure is preserved. At temperatures T ^ m/27i the spatial asymptotics of the 
Euclidean correlator reflect solely on screening. The tails no longer encode properly the 
real time singularity structure. This is more dramatic in the temporal direction, where 
the free fermion behaviour is reached at high temperature, even in the presence of a mass 
gap. 

Besides the asymptotics of the Euclidean correlator at finite temperature, one can 
probe its zero frequency (static) or zero momentum part, as first suggested by DeTar EI 



Indeed, from ( |6.4D we see that the large x^-behavior of the static part of the Euclidean 
correlator. 



= 0,x') = - f±;^,..'r'-_p(-yi (6.10) 

J 2,71 UJ 

reflects on the possible singularities of the spectral function in the pseudoscalar channel 
and/or screening. In the Schwinger model, the covariant dispersion relation u"^ = p^'^ + m'^ 
holds even at finite temperature, tying the singularity in to the singularity in u;. In 
general, this is not true and we would expect the singularities (poles, resonances and 
thresholds) to disperse in matter. The correlator ( p.lO| ) has been the object of interesting 
lattice simulations in the context of QCD [Q. Note that the zero momentum part of 
( |6.10| ) relates directly to the pseudoscalar susceptibility 



oo 



Xp = PEiu;„ =0,p'=0) = -2 I duj^^^, (6.11) 

Jo UJ 

where we have used the fact that (Jp{uj, 0) is and odd function of u at zero momentum]^. 
Susceptibilities of the type (|6.11|) reflect on the bulk properties of the theory at finite 



temperature . They are important dynamical parameters in the theory of linear 

response. For completeness, we remark that the zero momentum part of the Euclidean 



^The zero temperature limit ( 2.11 ) should be multiplied by sgn((7") to allow for antiparticles, and is 
consistent with this property. 
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correlator at finite temperature is given by 



f°° duj / ,a ^^ o-pfco" 0) cosh (a;(r — 1 1 



27r V ^ l-e^^ Jo sinh(^) 

(6.12) 

For high temperatures, the I.h.s. is dominated by the free field behaviour which should 
reflect itself in ap{uj ~ T,0). One can think of a numerical inversion of ( |6.12| ) (inverse 
Laplace transforms) as a way to probe the time-like character of (Jp{uj, 0) at zero mo- 
mentum. In the Schwinger model, the validity of the inversion can be checked term by 
term. This inversion is worth pursuing on the lattice for the finite temperature Euclidean 
correlators in QCD. 



7. Conclusions 

We have explicitly constructed the zero and finite temperature correlation functions 
for the scalar and pseudoscalar channels in the context of the Schwinger model. The finite 
temperature spectral functions have been obtained directly from the retarded correlators 
using operator methods, and from the Euclidean correlators by analytical continuation. 
The two procedures are in agreement. 

The temperature dependent spectral function in the Schwinger model reflects on ther- 
mahzed, but free mesons of mass m = e/y^. The poles and thresholds are temperature 
independent, with temperature dependent strengths. The absence of interactions in the 
bosonized version, makes the above results specific. In four dimensions and for theories 
such as QCD, we would expect the singularities to depend explicitly on the temperature 
(pions interact in the heat bath). However, we do not have a proof. 

Our expressions for the temperature dependent spectral functions can be written as 
the sum of all possible cuttings of free scalar diagrams with an arbitrary exchange of 
pseudoscalar mesons, for both the scalar and pseudoscalar channel. In their phase space 
version they are directly amenable to four dimensions. They provide compact expressions 



for cutting rules in real time [|IJ] in the context of theories. 
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The exact expressions for the Euchdean correlators indicate that the screening lengths 
in the Schwinger model are amenable to the meson mass m = e/i/vr and independent 
of the temperature, whatever the temperature. This is so, even though the meson is a 
composite of two fermions. We strongly suspect that this result is due to the point like 
character of the meson in the massless Schwinger model, and not just the anomaly. Indeed, 
we conjecture that in the massive Schwinger model, the screening lengths asymptote 2'kT 
at high temperature. We also conjecture that the heavy-light systems display screening 



lengths that asymptote vrT in this model [2S 



We note that in the Schwinger model the screening mass and the dynamical mass 
(leading singularity in the pseudoscalar channel at zero temperature) are temperature 
independent. The reason is that the process of mass generation in the model is through 
the Ua{^) anomaly, thus protected from matter corrections at all temperatures. This is 
not the case for physical excitations in four dimensions, to the exception of the t]' mass. 
This raises an interesting thought for four dimensional QCD. Could it be that the rj' mass 
survives matter effects? The standard scenario based on instanton-driven anomaly would 
not support this, as the instantons become suppressed at high temperature because of 
Debye screening. We remark, however, that the 1]' is not point like in nature. This point 
and the analogy are worth investigating, both in instanton and non-instanton approaches 
to the 7]' problem in matter. 

The above considerations suggest that possible criteria for a weakly interacting gas 
of massless fermions at high temperature might be formulated in terms of the sums, 
S = Se{x) + Pe{x), or the ratios, 71 = Se{x)/ Pe{x), taken along the spatial or temporal 
direction. For free massless fermions: S = and TZ = —1. These criteria hold in the 
presence of isospin, for a variety of combinations with proper isospin weighting. They 
might still be evaded by the Coulomb bound state approach following from dimensional 
reduction as discussed by two of us @]. This point is worth clarifying. 

The massless Schwinger model through its severe infrared sensitivity, provides a coun- 
terexample to the general lore on the spatial screening lengths based on dimensional 
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reduction in three and four dimensions for both QED and QCD. Ahhough the results are 
specific to the model, we think that they are still interesting. They provide a particular 
realization of field theoretical ideas at finite temperature, and show how naive arguments 
could be upset by infrared singularities. 

In four dimensions, QCD is also expected to suffer from strong infrared singularities 
at finite temperature. However, lattice simulations |l| suggest that at high temperature 



the screening lengths are compatible with the dimensional reduction scheme [Q, to 
the exception of the pion and its chiral partner. To what extent these results reflect on 
a free spectral function is not clear though, since these results are still compatible with 



correlations in the heat bath [^, ^]. 

Our arguments in the Schwinger model indicate that at low temperature T -C m/2iT, 
the large spatial asymptotics of the Euclidean correlators determines the leading isolated 
singularities of the spectral function, provided that the pre-exponents are properly iden- 
tified. This identification is rapidly lost as the temperature is increased, poles and cuts 
are hardly discriminated for T ~ m/2iT. If the analogy with the Schwinger model holds, 
the break up temperature in the QCD context would be T ~ 111^^/271 ~ 20 MeV. A 
similar remark holds for the behaviour of the Euclidean correlators along the temporal 
direction, except for the fact that the latter asymptote the free fermion behaviour at high 
temperature even in the presence of a mass gap. This point suggests that the lattice 
calculations of temporal correlators at high temperature may be misleading on the 
interacting aspects of the theory, since they reflect primarily on the free field content of 
the correlators. 

An interesting way to probe the dynamical part of the spectral function at zero mo- 
mentum is to try a numerical inversion of ( |6.12D for QCD. Also, if the lattice results could 
be analyzed in frequency space, then perhaps one may think about analytical continu- 
ations a la Baym-Mermin [^, assuming that one can reconstruct singularities from a 
numerical output on a discrete mesh ujn = 2TmT. In QCD at high temperature, it may 
be that the singularity is a simple cut throughout g° (analyticity aside) following the 
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change to an ionized (plasma) state. This simple fact is still elusive, however. Progress 
in these directions may come from real time alternatives to the present Euclidean time 
simulations. It should be emphasized, however, that the Euclidean correlators at finite 
temperature, are interesting in their own right. They reflect on spacelike physics in the 



form of screening and susceptibilities |24 



Finally, it would be interesting to see, whether finite temperature simulations of the 
massless Schwinger model, reproduce the exact results derived in this paper. Also, it 
would be amusing to test the ideas of finite temperature QCD sum rules |27| and finite 
temperature instanton (caloron) calculations |^ in this exactly solvable model. The 
extension of this work to higher correlation functions will be discussed elsewhere [ES . 



Note Added 

After finishing up this work we became aware of a recent paper by A. Smilga "Instan- 
tons in the Schwinger Model", NSF-ITP-93-151, where the scalar correlator is evaluated 
at finite temperature using instanton techniques. 
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Appendix A: The Euclidean fermionic Green's function 



In this appendix we calculate the fermionic Green's function defined by 

ipS{x,y) = t6{x-y)-P{x,y), (A.l) 

where ^is the Dirac operator in the presence of a U{1) gauge field 

p=^,,{d,,-ieA^), (A.2) 

and P is the projection operator onto the space of zero modes. We are interested in the 
case when there are no zero modes (P = 0). Since there are always zero modes present 
whenever the fiux $ 7^ we set $ = 

A = — ^ + ^10, 

A, = ^-9o0, (A3) 

where h^, hi are defined modulo one. The Dirac operator has two zero modes (one of 
each chirality) at (/io,/ii) = (1/2,0). However, these zero modes are cancelled by the 
determinant and do not give rise to singularities in the integration over ho and hi. In this 
section we will stay away from this point. Notice that if we write 

S {x, y) = e-T5*(^)f/ (x) g {x, y) (y) e-^5<A(j/)^ (A.4) 

where 

/ N f x^hn x^h^ \ / , X 

U (x) = expz27r f + j , (A. 5) 



then it follows from ( [A.l|) that 

i^g{x,0) =i5{x). (A.6) 



Because of translational invariance of the Green's function we have set y = 0. The 
boundary conditions on the Green's functions follow from the boundary conditions on the 
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fundamental fields^. Specifically 

g{x'> + P,x';0) = -e-^2.ho^ (a;; 0) , 

^(x°,x^ + L;0) = e-'^^^'g{x;0) . (A.7) 
From 75S'75 = —S it follows that S±^± = g±^± = 0. Equation ( |A.6|) then reduces to 

-2f3d.g+^{z) = (A.8) 

where z = {x^ + ix^) /(3. The above equations tell us that g-+{9+-) is a holomorphic 
(anti-holomorphic) function of z with a single pole at z = with residue l/27r/5 (— l/27r/?). 
The solution that satisfies the boundary conditions (|A.7|) can be written down at once 



9-+iz) = 7^—^ ] , = TT-^^^^'^iho + -,hi, z, 1, r), (A.9) 

2mp ^ z + m + riT 2mp 2 

where the function Ser was introduced and studied extensively by Kronecker. It can be 
expressed in terms of theta functions resulting in 



^^^^'^ - 27^^(3 Uz)Mhi - ^{h, + \)r) ' " ^^'^^^ 

The function g+-{z) is just the complex conjugate of The complete Green's function 
follows immediately by combining this result with the gauge transformation 

S^+{x,y) = g^+{z)e-^'^^^+^'^^h^^+^^ . (A.ll) 

The function g^^{z) can also be constructed directly in terms of theta functions by writing 



down a ratio with the boundary conditions ( |A.7D and the pole structure as dictated by 



( |A.10| ). The generalization of (|A.11|) to a nonvanishing flux sector can be found in . 



^Recall that $ = for the gauge fields, so that the gauge transformation one picks up under translation 
in the space direction is 1. 
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Appendix B: Contribution from the k = 2 sector. 



From (|4.1| ), the correlation function C++(a;) in the k = 2 sector reads 



:b.ii 



where ipi and 1IJ2 are the zero modes ([4.4|) for p = and p = 1, respectively. To evaluate 
first, we re-express ipi aiid 1^2 in terms of the 6*00 and 6*10 functions, respectively, 
using the following relations 



e 



e 

a + 



a 
b 

1 



(0,2ir) = exp(-27rra^ + 27rm6)6'oo(& + 2mr, 2ir), 
(0, 2iT) = exp(-27rra^ + 2mab)6io{b + 2mr, 2ir). 



:b.2) 



We are using (10) for (|0), (11) for (|^), etc. Second, we apply the Caspari-Kronecker 



formula 31 



eoo{2v, 2ir)9io{2u, 2ir) - ^10(21;, 2tT)9oo{2u, 2it) = 9n{v + u, iT)9u{v - u, ir) (B.3) 

Note, that in our case the argument {v — u) is /i- independent, so the remaining integrand 
\9ii{u + f,-jr)p over the harmonic fields is Gaussian. Thus 



C4 



1 1 



9iiiz, zr)Pe < e 



20(x)+2fli(O) 



>S 



X 



(B.4) 



where t = xq/ P + hi + ir^xi/L + 1/2 — ho). Since the integration over hi reduces the 
double sum in \9ii{t,iT)\'^ to a single one, the correlator reads 



e-'-^\9ii{z,tT)\'e-V^ < e2<^(-)+20(O) >^ y ^^^^-2nr(n+l-ho + ^) 

^ Jo 



1 1 ■^^^ 



C++(x) 

where z = (xq + ixi)/ (3 denotes the relative coordinate 



^ ^ ^e-^|^n(^,^r)|2e-|^ < e^^^^^+^^W >5, (B.5) 
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Figure Captions 

Fig. la. Temperature effects on tfie sources (a). Tlie dasfies are temperature insertions. 

Fig. lb. Diagrammatic expansion of the retarded pseudoscalar correlator Pe- The 
straight hnes are retarded propagators A^^, and the broken hnes the imaginary part of 
the Feynman propagator ImAp. 

Fig. Ic. The same as (b) but for the retarded scalar correlator Se- 

Fig. 2a. Leading processes contributing to the scalar spectral function ug in the heat 
bath. 

Fig. 2b. Leading processes contributing to the pseudoscalar spectral function up in the 
heat bath. 
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